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We discuss properties of entanglement measures called I-concurrence and tangle. 
For a bipartite pure state, I-concurrence and tangle are simply related to the purity 
^| ■ of the marginal density operators. The I-concurrence (tangle) of a bipartite mixed 

. state is the minimum average I-concurrence (tangle) of ensemble decompositions of 

, pure states of the joint density operator. Terhal and Vollbrecht [Phys. Rev. Lett. 

85, 2625 (2000)] have given an explicit formula for the entanglement of formation of 
isotropic states in arbitrary dimensions. We use their formalism to derive comparable 
expressions for the I-concurrence and tangle of isotropic states. 
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> 

Entanglement among quantum systems is a distinctive feature of quantum mechanics |y, 
H, || and an indispensable ingredient in various kinds of quantum information processing 



O 

oo 



protocols [|, U, e.g., quantum teleportation. It is desirable to have a general theory of 
entanglement, but though important results have been obtained, a general theory has proven 
elusive because of the complex nature of entanglement for all but the very simplest joint 
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quantum systems. Several measures of entanglement have been suggested and investigated: 
the entanglement of formation || [7|, the entanglement of distillation |6|, [7|], the relative 
entropy || |9|, the robustness of entanglement [[KJ], and others. In this paper we explore 



Wootters [11, 12 1 for pairs of qubits and generalized to arbitrary bipartite systems by Rungta 



a measure of entanglement called the I-concurrence, which was introduced by Hill and 



et al. |TJ|. Before turning to the I-concurrence, we summarize here relevant results from 



entanglement theory which will facilitate the discussion of I-concurrence. 



Work on entanglement has served to identify certain a priori axioms for a good measure 



of entanglement [TP1. Entanglement characterizes intrinsically quantum- mechanical correla- 



tions between quantum systems. If E(p) denotes an amount of entanglement in the joint 
state p of several quantum systems, a fundamental requirement for E to be a good measure 
of entanglement is the following: 

1. E(p) does not increase, on average, under local operations and classical communication 



(LOCC). Such a measure is called an entanglement monotone |L4 



An entanglement monotone must remain the same under the action of reversible LOCC. 
Hence, it is invariant under local unitary transformations since these can be locally reversed. 
Moreover, any state p can be converted into any separable state using LOCC; therefore, an 
entanglement monotone takes a common minimum value for all separable states, which 



can always be adjusted to zero [14]. Thus we can impose an additional, nonnegativity 



requirement for a good entanglement measure: 

2. E(p) > and goes to zero if and only if p is a separable state. 
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Henceforth, we restrict our discussion to bipartite quantum systems. In three or more 
Hilbert-space dimensions, more than one measure is required to characterize the entangle- 
ment of bipartite states [ID, 11, I5| : to specify completely a joint pure state of a d x d system, 
up to local actions, requires d — 1 independent Schmidt coefficients. Vidal [[14]] showed how 
to construct an infinite family of entanglement monotones for bipartite systems. This family 
of entanglement measures, denoted by ^(p), is characterized by the following conditions. 

(i) For a pure state | | , the measure is a function of the marginal density operator pa = 
tr#(|\I/) (\I/|), i.e., = /(pa), where the function / has the following properties: 
(a) invariance under unitary transformations, i.e., f(UpAU^) = /(pa)> which implies 
that /(pa) is a function of the eigenvalues of pa, and (b) concavity, i.e., f(\\Pi + 
A2P2) > Ai/(pi) + A 2 /(p 2 ), where Ai, A 2 > and Ai + A 2 = 1. 

(ii) For a mixed state p, the measure is defined to be the convex-roof extension of the 
pure-state measure, i.e., the minimum average value of the measure over all ensemble 
decompositions of p: 



u(p) = min 



Vidal []nj showed that if a measure belongs to the above family, it satisfies property 1 and 
thus is an entanglement monotone. The I-concurrence, as we show in Sec. ^|, belongs to this 
family. Vidal also showed that when restricted to pure states, any entanglement monotone 
satisfies condition (i). 

A privileged example of an entanglement monotone comes from choosing the unitarily 
invariant concave function of condition (i) to be the von Neumann entropy, /(pa) = S(pa) = 
— tr (pa log 2 Pa)- This entanglement measure plays a special role because of the asymptotic 
reversibility of the processes of entanglement dilution and concentration for pure states, 
with von Neumann entropy entering as the currency for these reversible transformations 
16fl : nS(pA) Bell states can be converted by LOCC to n copies of in the limit as n 
goes to infinity and vice versa ||. The corresponding entanglement monotone for pure 
states, Ef(^) = S(pa), when extended to mixed states by the convex roof, is called the 
entanglement of formation || [/]]: 



E f (p)= min \y^p j E f (^ 1 



Unlike the situation for pure states, however, it is not known whether the entanglement of 
formation gives the asymptotic cost in Bell states for creating many copies of a mixed state p 
using LOCC. This identification depends on whether Ef(p) is additive 0, 0. Though Ef(^) 



is additive on pure states, this only implies that Ef(p) is subadditive on mixed states. If Ef(p) 
is strictly subadditive, then some other entanglement monotone, given by the von Neumann 
entropy for pure states, but less than Ef(p) for mixed states, quantifies the asymptotic cost 
of creating many copies of a mixed state @, [17] . 

The special role of von Neumann entropy has been reconciled with the existence of many 



other entanglement monotones by Vidal [|14[] and Nielsen |18 |. Vidal considered a class 



of entanglement monotones that are additive on pure states and showed that they have 
different asymptotic properties than the entanglement of formation. Nielsen formulated a 
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precise property of asymptotic continuity that is satisfied by the von Neumann entropy and 
showed that this property is the key to relating entanglement monotones to the asymptotic 
reversibility of entanglement dilution and concentration that holds for pure states. Indeed, 
Nielsen showed that any entanglement monotone that is additive on pure states and that 
satisfies his asymptotic continuity property is given, for pure states, by the von Neumann 
entropy of the marginal density operator. This result of Nielsen's leaves open the question 
of the entanglement cost for mixed states. 

Entanglement monotones other than the entanglement of formation, though they are not 
related to asymptotic transformations, are nonetheless important for characterizing LOCC 
transformations between finite numbers of copies of states. Indeed, they have been used to 
quantify the probability and fidelity of exact and approximate LOCC transformations be- 
tween pure states [|l^, pi], ^TJ . In this paper we investigate the properties of entanglement 
monotones related to a particular entanglement measure called the I- concurrence. 

Hill and Wootters fill introduced the concurrence as a measure of entanglement for pairs 
of qubits. The concurrence for a pure state of two qubits is defined with the help of the 
qubit spin-flip operator, and it is extended to mixed states as the convex roof. Wootters 
@, [12| went on to derive an explicit formula for concurrence of an arbitrary joint state p of 
two qubits in terms of the eigenvalues of y/ppy/p, where p is obtained by spin flipping p, and 
showed how to calculate the corresponding entanglement of formation from the concurrence. 



Rungta et al. [|T^] generalized the notion of concurrence to pairs of quantum systems of 
arbitrary dimension. This generalized concurrence for a joint pure state of a oIa x &b 
system is simply related to the purity of the marginal density operators: 



C(¥) = p[\ - tr(pi)] = p[l - tr(p|)] . (1.3) 
The generalized concurrence is known as the I-concurrence because it is defined in terms of 



the universal- inverter superoperator [[13], which was shown to be a natural generalization to 
higher dimensions of the spin flip for qubits. For the purposes of this paper, we need not 
introduce the universal inverter, simply taking Eq. (|1.3|) as the definition of the I-concurrence 



for pure states. The I-concurrence is extended to mixed states p by the convex roof: 

£p^><^| = 4 . (1.4) 



C(p) = min iy"p,C(^ 1 



For pairs of qubits, Coffman, Kundu, and Wooters |22| called C 2 (p) the tangle of the 
state p. In this paper we prefer to refer to C 2 (p) simply as the squared I-concurrence and 
to reserve the term tangle for the quantity obtained by extending t(\&) = C 2 (\l/) to mixed 
states by the convex roof: 



t( P ) = min \ J2PjC 2 (^j 



Osborne |23[ has calculated the tangle of rank-2 density operators in 2 x d systems. As noted 
by Osborne, the definition of tangle in Eq. ( |1.5|) does no violence to the original usage of 
Coffman, Kundu, and Wootters, because r(p) = C 2 (p) for bipartite qubit states, the reason 
being that all pure states in the optimal ensemble decomposition have the same concurrence. 
Osborne found the tangle as defined in Eq. ( |1.5|) to be a more natural quantity in his analysis 
of rank-2 density operators, just as we find it to be more natural for isotropic states in our 
analysis below. 
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The motivation for this paper is to show that, like the entanglement of formation, the 
I-concurrence and tangle are entanglement monotones (Sec. |2|) and then to provide explicit 
formulae for the I-concurrence and tangle of isotropic states (Sec. |3|). For our analysis of 
isotropic states, we use the formalism developed by Terhal and Vollbrecht [23J] to calculate 
the entanglement of formation for isotropic states. We find that the tangle for isotropic 
states is closely related to the corresponding entanglement of formation (Sec. H). 



2. PROPERTIES OF I-CONCURRENCE AND TANGLE 

Henceforth, we will omit the "I" when referring to the I-concurrence. Consider the 
Schmidt decomposition of an arbitrary joint pure state |^} of a x ds system: 

l^> = E y/ih\ a i) ® \ h i) = Ev^W ® U B\ e j) , (2- 1 ) 
j j 

The squared Schmidt coefficients, pj, are the eigenvalues of the marginal density operators, 
Pa and ps, of the two systems, and the vectors |a_j) and \bj) make up the orthonormal bases 
that diagonalize the marginal density operators. These bases are connected to a fiducial 
orthonormal basis {|e-,)} by unitary transformations Ua and Ub- The state |\&) is specified 
by its Schmidt vector /x = (y/JIi, ■ ■ ■ , y/Pd) and the unitary operators Ua and Ub- 

The tangle, or squared concurrence, of the pure state |^) [see Eq. (|1.3| )1 is given in terms 
of the Schmidt coefficients by 

r(¥) = C 2 W = 2 (l - £ A = 4 ^ Wk = CV) • (2.2) 

C 2 (\l/) is conserved under by local unitary transformations because it is a function only of 
the Schmidt coefficients [property (i)a]. It varies smoothly from 0, for pure product states, 
to 2(d — l)/d, where d = min(d^, ds), for maximally entangled pure states. 

That C 2 (\l/) is a concave function of pa (pb) [property (i)b] follows from the fact that 
f(x) = —x 2 , < x < 1, is a concave function, since for any concave function /, the mapping 
p — > tr(/(p)) is concave [ffiSfl). To see this, let p = X1P1+X2P2, with Ai, A2 > and A1+A2 = 1, 
and let p = J2j Pj\4>j) (4>j\, Pi = Efe9fe|6)(Cfc|: and p 2 = J2i ri\r]i) (m\ be eigendecompositions. 
Then we have 

tr(/(p)) = 

3 

= I]/( A l(0ilPll0i) + ^(<f>j\P2\(f>j)) 
j 

> E A i/«^IPil^» + W«^|pa|^» 

3 

= Ai E / (E I (0iie*> I a ft) + A 2 E / (E I (htotfri) 
j v fc j v 1 

> Ai E I to l&> iVfe) + A 2 E I (fah) \ 2 f(ri) 

j,k 3,1 

= Ai £/(«?*) + AaE/(n) 

= A!tr(/( Pl )) + A 2 tr(/(p 2 )) , (2.3) 
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where the two inequalities follow from the concavity of /. Since adding a constant and 
multiplying by a positive constant doesn't change concavity, C 2 (^>) is a concave function 
of pa- Since C 2 (^/) satisfies property (i), its extension (|1.5| ) by the convex roof to give the 
tangle r(p) of mixed states is an entanglement monotone. 

The square root being an increasing concave function, it preserves concavity. Thus the 
concurrence C(^) is also a concave function of pa, and its extension (|1.4 ) to give the 
concurrence C(p) of mixed states is an entanglement monotone. Notice that since C(\l/) is 
zero if and only if |\&) is a pure product state, r(p) and C(p) are zero if and only if p is 
separable. 

From the properties of an entanglement monotone or directly from the convex-roof con- 
struction, it is clear that the concurrence C(p) is a convex function of bipartite density 
operators. Since the square is an increasing convex function, it preserves convexity, thus 
making the squared concurrence, C 2 (p), a convex function of bipartite density operators. 
As a convex function that agrees with the convex roof r(p) on pure states, the squared 
concurrence is guaranteed to satisfy C 2 (p) < r(p). We are unable to say whether C 2 (p) is 
itself an entanglement monotone, although it seems unlikely that it satisfies the property of 
not increasing under local measurements. 



3. TANGLE AND CONCURRENCE OF ISOTROPIC STATES 

3.1. Isotropic states 

In this section we derive the tangle and concurrence of isotropic states. At the end of 
the section, we compare our results with the those for the entanglement of formation of 
isotropic states obtained by Terhal and Vollbrecht []24|] . The tangle of isotropic states shares 
important features with the entanglement of formation, but the concurrence is significantly 
different. 

Isotropic states are a class of mixed states for d x d systems (two qudits); they are convex 
mixtures of the maximally mixed state, = I®I /d 2 , with a maximally entangled state, 

l*+>= * 5>,.>®| ej .>. (s- 1 ) 

V a j=1 

Such mixtures can be expressed as 

Pf = (I - l* + X* + l) + F|*+}<*+| , (3.2) 

where F = (^/ + \pf\^ + ) , satisfying < F < 1, is the fidelity of p F and |^/ + ). These states 
were shown to be separable for F < 1/d |2B|, ^7 1 . 



The isotropic states are special in the sense that they are invariant under the action of 
the twirling superoperator T 



T(p F ) = J dUU® U*p F U ] ® U* j = p F . (3.3) 
Indeed, the twirling superoperator reduces any two-qudit state p to an isotropic state 



^(P) = PF(p) , 



(3.4) 
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where F(p) = is the fidelity of p and |^ + ). Twirling the pure state of 

Eq. ( |2.1|) yields an isotropic state 



T(|*)(*|) = Pi?(M! 



v) 



where the fidelity is given by 



X \f^k{ e i\UA\ek) (ej\U B \e k ) 

2 

xi \fpk{UA)jk{UB,)jk 

2 



kk 



V = UbU\ being a unitary matrix. It is easy to see that 

F^,V)<F(fi,I), 
since \V k k\ < 1; equality holds if and only if V = I exp(i5). 



(3.5) 



(3.6) 



(3.7) 



3.2. Reduction to a single extremization 

We turn now to finding the tangle and concurrence of isotropic states, using the technique 



developed for entanglement of formation by Terhal and Vollbrecht [24]. We proceed with 
the analysis in terms of the tangle. The reader should note that nothing would change in 
the analysis if we replaced pure-state tangle C 2 (\I/) with pure-state concurrence C(^), thus 
analyzing C(p) instead of r(p), till we get to Sec. |3]^. At that point, we explicitly note the 
differences between the results for tangle and concurrence. 

The tangle of an arbitrary bipartite state p satisfies the following inequality, 



r(p) 



= min {X^C 2 ^) 



X^X^I 



p 



^p j F( H ,V j ) = F(p) 



(3.8) 



(3.9) 



The inequality in Eq. ( |3.9| ) follows because an optimal decomposition of p, i.e., one which 
achieves the minimum in Eq. ( |3.8| ), automatically generates a set {pj, Hj, Vj} that satisfies the 
constraint in Eq. (|3.9| ). In contrast, a set {pj, fj,j, Vj} that achieves the minimum in Eq. ( |3.9|) 
generally does not generate an ensemble decomposition of p. Additional simplification is 
achieved by splitting the minimization in Eq. (|3.9j) into two parts, i.e., 



r(p) > mm .fepjCfyiFj) ^ptf = F(p)\ = r(F(p)) , (3.10) 
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with 

C 2 (F)=miJc 2 (») F(/i,V) = ± 



(3.11) 



The function r(F) defined in Eq. ( |3.10| ) is a function of the single parameter F = F(p). 

We can reduce the minimization problem further by noting that if v is the vector of 
Schmidt coefficients that provides the minimum for C V (F), then F' = F(u, I) > F(u, V) = 
F [Eq. flO])] and 

C 2 {F') = Cj(F') < C 2 (y) = C 2 V (F) . (3.12) 

We find an explicit expression for C 2 (F) below and show that it is monotonically increasing, 
from which it follows that 

C V (F) > C 2 {F') > C 2 {F) . (3.13) 
Applying Eq. ( &T3j ) to Eqs. ( PH) and (|3~TT|) yields 

t(F)= mmj^^ 2 ^) /]• (3-14) 

C 2 (F)=min|c 2 ( M ) ^(Ev^) 2 = ^ (3-15) 

Notice that the inequality C 2 (F') > C 2 (F) in Eq. (|3.14|) requires us to assume that F > 1/d, 
since the minimum that defines C 2 (F) does not exist for F < 1/d. This is not a problem 
for the analysis of isotropic states, since isotropic states with F < 1/d are separable, having 
t(pf) = and C(pf) = 0. The function t(F) is by definition a convex function of F. 
Indeed, t(F) can be defined as the largest convex function that is bounded above by C 2 (F); 
it is given either by C 2 (F) or by straight-line segments that connect points on the graph of 
C 2 {F) and lie beneath C 2 {F). 

We have shown that the tangle for any bipartite pure state is bounded below by r(p) > 
r(F(p)). Now we use the twirling superoperator to show that isotropic states achieve this 
lower bound. To do so, consider an isotropic state pp- Let {pj, Fj, /-ij} be a set that 
achieves the minimum in Eqs. (|3.14|) and (|3.15 ) with this value of F, and let {|^}} be 



states constructed from the Schmidt vectors {^} with Vj = I. The state formed from this 
ensemble, 

P = EPil*iX*il. (3-16) 

j 

has a tangle that satisfies r(p) < J2jPjC 2 (^j) = t(F). Twirling p gives 

T(p) = Y^p/K M/ a (vi, ; j = J2 Pj p F . = p F , (3.17) 

j j 

where we have made use of Eq. (|3T5|) to write T(\tyj)(tyj\) = p F ^. Since the local operations 
involved in twirling cannot increase the tangle, we have r(p) > r(p^), from which follows 
the upper bound r(p F ) < t(F). 

Combined with the lower bound r(p F ) > t(F), this shows that the tangle of an isotropic 
state p F is given by the function t(F). We have thus reduced the problem of finding the 
tangle of an isotropic state to a single minimization, that of finding the function C 2 (F), from 
which t(pf) = t(F) can be constructed as described above. If we follow through the steps 
of this subsection for the concurrence, instead of the tangle, we find that the concurrence of 
an isotropic state, C(p F ), is the largest convex function that is bounded above by C(F). 
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3.3. The extremization 



Following the method of Terhal and Vollbrecht p4[ , we calculate C 2 (F) using the method 
of Lagrange multipliers; i.e., we minimize C 2 (/x) [Eq. ( [2.2|) 1 subject to the constraints 



j 



'Fd 



(3.18) 
(3.19) 



with Fd > 1. In doing the extremization, we allow for the possibility that the minimum 
might not have all nonzero Schmidt coefficients by explicitly considering all the cases where 
the number of nonzero Schmidt coefficients varies from 1 to d. The condition for an extremum 



is 



(3.20) 



where Ai and A2 are Lagrange multipliers. The three solutions of this cubic equation for 
y/JIj sum to zero, so there are at most two real, positive solutions. Letting 7 and 5 denote 
these two positive solutions, with 7 > S, the possible Schmidt vectors fx have coefficients 




l,...,n, 

n + 1, . . . , n + m, 
n + m + 1, . . . , d, 



(3.21) 



where n + m < d and n > 1. The corresponding extrema of C 2 (fj,) are 



C 2 nm (F) = 2 (l - n 7 4 - m5 4 ) 



(3.22) 



with the constraints 



ny 2 + m5 2 
wy + m5 



'Fd 



Solving Eqs. ( p.23|) , we obtain 



Trim 

% m (F) 



n\/ Fd ± J nm(n + m — Fd) 
n(n + m) 

V^Fd - n>y± m mVFd =f Jnm(n + m- Fd) 



m 



m(n + m) 



(3.23) 

(3.24) 
(3.25) 



The relation 5^ m (F) = JmniF) means that as we vary n and m over all possible values, we 
need only consider the upper sign; henceforth we drop the signs, using always the upper 
sign. For the expressions ( |3.24j ) and (|3.25|) to give real values, the quantity inside the square 
root must be nonnegative, which implies that Fd < n + m; moreover, in order that 5 nm be 
nonnegative, we must have Fd > n. It is easy to see that 8 nm (F) < V Fd/ (n + m) < y nm (F), 
confirming that the choice of the upper sign corresponds to our assumption that 7 > 5. 
Notice also that n = is ill defined, as expected, thus requiring n > 1. 

The function C 2 (F) we seek is the minimum of C 2 m (F) over all choices of n and m. 
We can perform the minimization explicitly by regarding n and m as continuous variables 
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(here our treatment departs from that of Terhal and Vollbrecht). The task is to minimize 
C 2 m (F) on the parallelogram defined by 1 < n < Fd and Fd < n + m < d. Notice that the 
parallelogram collapses to a line when Fd = 1, i.e., at the separability boundary. As already 
noted, within the parallelogram we have 7„ m (-F) > S nm (F) > 0; ^ nm (F) = 5 nm (F) if and 
only if n + m = Fd, and S nm (F) = if and only if n = Fd. We first calculate the derivatives 
of 'jnm(F) and S nm (F) with respect to n and m by differentiating the constraints ( |3.23| ): 

c?7 1 2^5 — 7 2 
dn 2n 7 — 5 ' 
05 1 7 2 



dn 2m 7 — 5 

05 1 2 1 5-5 2 

dm 2m 7 — 5 ' 
9 7 1 5 2 
dm 2n 7 — 5 

These can be used in Eq. ( |3.22| ) to calculate the partial derivatives of C 2 m (F): 



(3.26) 



dC 2 

dn 
dC 2 



2 7 2 [7 2 -25( 7 + 5)], (3.27) 

25 2 [5 2 -27(7 + 5)] < -6<5 4 < . (3.28) 
dm 

It is useful to introduce coordinates u = m — n and v = m + n, which correspond to motion 
parallel to and perpendicular to the m + n = constant boundaries of the parallelogram. The 
derivative of C 2 with respect to u is 

dC 2 1 

^- = --( 7 + 5)(7-5) 3 <0. (3.29) 

The inequalities in Eqs. ( |3.28| ) and ( |3.29| ) follow immediately from the fact that 7 > 5 > 
within the parallelogram. It should be clear that dC 2 m /dm is strictly negative within the 
parallelogram and dC 2 m /du is strictly negative except on the boundary m + n = Fd, where 
it is zero. These results mean that the minimum of C 2 m (F) occurs at the vertex n = 1, 
m = d — 1, thus giving 



C\F) = Cl^iF) = 2(l - 7^-1 - (d - lKVj • ( 3 - 3 °) 



Here 



7M _i(F) = JF/d (l+wVd^l) , (3.31) 



<f M _i(F) = ^F/d(l-w/Vd^T) , (3.32) 



with u> = y(l — F)/F. Henceforth we omit the subscripts that specify the case n = 1 and 
m = d — 1, this being the only case of interest. 

We need to confirm that C 2 (F) is monotonically increasing. Differentiating the con- 
straints ( |3.23| ) gives 
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from which we calculate 



dC 2 



4W-75(7 + ^) > 
b 



(3.35) 



equality holding if and only if 5 = 0, which requires Fd = 1, i.e., the separability boundary. 
We conclude that C 2 (F) [and C(F)] are monotonically increasing for 1/d < F < 1, as 
promised. It will be useful below to write dC 2 /dF explicitly in terms of F and d: 



dC 2 
dF 



8^ (l + wVd^V) (l 



w 



s/d 



l + w- I Vd 



1 



(3.36) 



3.4. Results 



3.4-1- Two qubits 

The case of two qubits (d = 2) is special, so we discuss it separately. For two qubits, 
there is one extremum, C 2 (F) = Cf^F) = (2F- l) 2 , 1/2 < F < 1. Since C 2 {F) is a convex 
function of F, it follows that the tangle of an isotropic state pp is given by 



t(pi 



0, 

(l-2Ff 



< F < 1/2, 
1/2 < F < 1, 



(3.37) 



and that the pure states in an optimal ensemble decomposition for the tangle all have the 
same tangle. The same conclusions hold for the concurrence, C{pp) = C(F) = 1 — 2F, 
1/2 < F < 1, which agrees with the general expression derived by Wootters [J2j . 




FIG. 1: Plots of Cl x (F) (dashed line), Cf^F) (dotted line), and C 2 2 (F) = C 2 (F) (solid line) for 
d = 3. 



5.^.^. Two qutrits 

We consider the case of two qutrits (ci = 3) separately as an illustration of what happens 
in the general qudit case. Of the three extrema, Cf^F), C 2 X (F), and C 2 2 (F), we already 
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FIG. 2: First (solid line) and second (dotted line) derivatives of C 2 (F) = C 2 2 {F) for d = 3. 

know that the minimum is given by Cf 2 (F) } i.e., C 2 (F) = Cf 2 (F), 1/3 < F < 1. This fact 
can also be seen from Fig. [T], where we have plotted the three extrema. 

To calculate r(p F ), it is necessary to analyze the behavior of C 2 (F). In Fig. ^ we 
plot the first and second derivatives of C 2 (F). The first derivative confirms that C 2 (F) 
is monotonically increasing. The second derivative changes sign from positive to negative; 
C 2 (F) changes from convex to concave where the second derivative vanishes. The tangle 
t(pf) is the largest convex function bounded above by C 2 (F), which is constructed in the 
following way. Find the (unique) line tangent to C 2 (F) that passes through the point 
(F = 1, C 2 = 4/3); for F smaller than the tangent point, the tangle is given by C 2 (F), but 
for F larger than the tangent point, the tangle is given by the line. The tangent point is 
found by solving the equation 8C 2 /dF = [4/3 - C 2 (F)]/(1 - F), which gives F — 8/9. The 
slope of the line is (dC 2 /dF)\p=8/9 = 3, and the tangle at the tangent point is C 2 (8/9) = 1. 
Thus the tangle for d = 3 is 

f , F < 1/3, 

t(pf) = { C 2 (F) , 1/3 < F < 8/9, (3.38) 

[3(F-l)+4/3, 8/9<F<l. 

This function is plotted in Fig. |3|. 

The behavior of C 2 (F) means that the pure states in an optimal ensemble decomposition 
for the tangle all have the same tangle for F < 8/9, but have two values of tangle, 1 and 
4/3 (maximal entanglement), for F > 8/9. 

Of the three extrema for the concurrence, Cn(F), C2i(F), and C\2(F), the minimum is 
given by Ci 2 (F) = C(F), a fact confirmed by the plots in Fig. [|. In contrast to the situation 
with the tangle, however, the second derivative of C(F) = Ci 2 (F) is always negative (see 
Fig. U), which means that C(F) is concave. Therefore the qutrit concurrence is the straight 

line passing through the points (F — 1/3, C — 0) and (F — 1, C — ^4/3), i.e., 

C{pF) = {°/3(F-l/3), 1/3<F<1. (3 - 39) 

The concavity of C(F) means that the pure states in an optimal ensemble decomposition 
for the concurrence are either product states or maximally entangled states. 
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FIG. 3: Tangle r{pp) for d = 3 (dotted), d = 10 (short-dashed), d = 100 (long-dashed), and 
d = 10 000 (solid). The solid line is indistinguishable from the asymptotic tangle. 




FIG. 4: Plots of C n (F) (dashed line), C 21 (F) (dotted line), and C X2 {F) = C{F) (solid line) for 
d = 3. 

3.4.3. Two qudits 

For arbitrary d, we have already established that C 2 (F) = C^ d _ 1 (F). It turns out that 
for arbitrary d > 3, C 2 (F) has the same behavior as for d = 3; i.e., it changes from convex 
to concave as F increases. We get at this behavior by calculating the second derivative of 
C 2 from Eq. Q3gg ): 

<9 2 C 2 1 d (dC 2 \ 6 d-2 / 2 d 2 -8d + 8 n 2 1 A 

- - = -— — — = ; 1 H w - 2w w . 

OF 2 2wF 2 dw\dFJ wd^d^l\ 3{d-2)y/d=T 3 J 

(3.40) 

Notice that for d = 2 this expression simplifies to d 2 C 2 /dF 2 = 8, as it should. For d > 3, the 
polynomial in large parentheses in Eq. ( p.40[ ) clearly has one root for positive w; it is easy to 
verify that this one root occurs in the range of interest, i.e., between F — 1/d (w — y/d — 1) 
and F — 1 (w — 0), showing that C 2 (F) changes from convex to concave as F increases. 
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— ' 1 1 1 1 1 ■- F 

0.4 0.5 0.6 0.7 0.8 0.9 1.0 



FIG. 5: First (solid line) and second (dotted line) derivatives of C(F) = C '12(F) for d = 3. 

We find the tangle r(pp) just as for d = 3, i.e., by finding the line tangent to C 2 (F) 
that passes through the point (F — 1, C 2 — 2(d — l)/d). Solving the equation dC 2 /dF = 
[2(d — l)/d — C 2 (F)]/(1 — F) gives a tangent point at F = 4(d - l)/d 2 , where the tangle 
and slope are C 2 = 2(2d - 3)/d(d - 1) and dC 2 /dF = 2d/{d - 1). Thus the tangle is given 
by 

r , f < i/d, 

t( Pf ) = I C 2 {F) , 1/d < F < 4{d - l)/d 2 , (3.41) 

[ 2d(F - l)/(d - 1) + 2(d -l)/d, 4(d - l)/d 2 < F < 1. 

The tangle is plotted for a few representative values of d in Fig. |3|. The pure states in an 
optimal ensemble decomposition for the tangle all have the same tangle for F < 4(d — l)/d 2 , 
but have two values of tangle, 2(2d — 3)/d(d — 1) and 2{d — l)/d (maximal entanglement), 
for F > 4(d — l)/d 2 . As d — > 00, r(pp) — > 2F becomes a linear function. 

The concurrence also has the same behavior generally as for d = 3. A calculation of 
d 2 C/dF 2 using Eqs. ( [3.36D and ( |3.40|) shows that C(F) is a concave function on the range 



of interest. Thus the concurrence of isotropic states is a linear function between zero con- 
currence at the separability boundary and maximal entanglement at F — 1: 

(0 F<l/d, 

C ^ F) = \pd/{d-l){F-l/d) , l/d<F<\. (3 - 42) 

The pure states in an optimal ensemble decomposition for the concurrence are either product 
states or maximally entangled states. As d — >• 00, C(pp) — > \f2F . 



4. CONCLUSION 

It is informative to conclude the paper by comparing our results for the tangle of isotropic 



states with the results of Terhal and Vollbrecht [24| for the entanglement of formation. 

If one follows through the procedure of Terhal and Vollbrecht J24j for finding the entan- 
glement of formation, one finds that the function C 2 (F) is replaced by the function 



£(F) = tf 2 ( 7 2 ) + (l-7 2 )log(d-l), 



(4.1) 
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where 7 = 7i,d-i and H 2 (x) = — xlogx— (1— x) log(l— x) is the binary entropy function. The 
entanglement of formation, E(pp), is given by the largest convex function that is bounded 
above by E(F). Terhal and Vollbrecht conjecture that E(F) has the properties that we 
find here for C 2 (F), i.e., that E(F) has a single inflection point in the range 1/d < F < 1, 
changing from convex to concave as F increases. Given this conjecture, the entanglement 
of formation is found by finding the straight line tangent to E(F) that passes through the 
point (F — 1, E — logd). Remarkably the tangent point for E(F) is the same as for C 2 (F), 
which gives the entanglement of formation as 

r , i/d < f, 

E f (p F ) = E(F) , l/d < F < 4(d - l)/d 2 , (4.2) 

[d\og(d-l)(F -l)/(d-2)+\ogd , A(d-l)/d 2 < F < 1. 

Asymptotically, as d — > 00, the entanglement of formation becomes Ef(p F ) — > Flogd. 
The similarity of the tangle and entanglement of formation must reflect a deep connection 
between the two, at least for isotropic states, but we have not been able to find a simple 
reason for this similarity. 
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